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ABSTRACT 

Using a model for the heat flux through a single 
drop, a computer simulation of the dropwise condensation 
process was carried out. For a given set of condensing 
conditions, the computer model predicts heat transfer 
coefficients and drop distributions which agree with 
experimental observations. Subsequently, an analysis was 
conducted with the model by varying the condensing para- 
meters of site density, saturation temperature, and tempera 
ture difference. From a comparison of the analysis with 
experimental results, the effects of the previously 
mentioned parameters were determined. 
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NOMENCLATURE 

a = constant defined in Eg. (A. 15) 

A = Area Cft^} 

b = constant defined in Eg. (A. 15) 

c = constant defined in Eg. (A. 15) 

d = constant defined in Eg. (A. 15) 

D = drop diameter (ft) 

h = heat transfer coefficient defined in Eg. (F.3) 
CBtu/hr ft2 «F) 

H = latent heat of vaporization (Btu/lbm) 

K = thermal conductivity (Btu/hr ft°F) 

M = molecular weight 

Q = rate of heat transfer (Btu/hr) 

r = radius of drop (ft) 

R = resistance (hr ®F/Btu) 

R = universal gas constant 

t = time or drop age (seconds) 

T = temperature 

V = specific volume (ft^/lbm) 

V = volvime of drop (ft^) 

a = mass transfer accommodation coefficient, 

condensation coefficient (dimensionless) 

AT = temperature difference (®F) 

3 

p = density (Ibm/ft ) 

a = surface tension (Ibf/ft) 
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NOMENCLATURE 

(Cont'd) 

SUBSCRIPTS 

c = curvature 

cyl = cylinder 

dc = drop conduction 

fg = refers to change by condensation 

g = vapor 

i = interfacial 

max - maximum 

min = minimum 

nuc = nucleation 

s = saturation 

t = total 
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CHAPTER I 
INTRODUCTION 

Over the past few decades condensers have remained 
essentially unchanged, while other components of the steam 
power plant have been improved. The potential for improvement 
of the condenser exists if the phenomenon of dropwise con- 
densation can be utilized. Dropwise condensation is the 
process where vapor, condensing on a cold and non-wetting 
surface, forms liquid drops which grow and coalesce and 
ultimately roll off. For a surface which the liquid wets, 
a continuous liquid film forms and covers the surface. 

This latter process is called filmwise condensation and is 
the basis of current condenser design. Experimentally, 
heat transfer coefficients for dropwise condensation have 
been found to be an order of magnitude larger than the 
coefficients for filmwise condensation [1]. With this 
superior heat transfer property, condensers can be made 
smaller by weight and volume. 

There are three problems retarding the use of 
dropwise condensation. First, there is today no well proven 
surface or surface coating which promotes the formation of 
drops. Secondly, the presence of noncondensable gases in 
the vapor reduces the heat transfer coefficient. There is a 
need for a means of reducing the noncondensables to an 
acceptably low limit. Finally, there is a lack of understand- 
ing of the effect of certain condensing parameters. This 
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study is to aid understanding by providing a theoretical 
prediction of the heat transfer properties of dropwise 
condensation by use of computer simulation of the growth and 
coalescence of drops from nucleation to the size at which 
the drops slide off the surface. 

This investigation models condensation on a 
vertical surface. Initially the surface is bare. Then 
steam, without noncondensable gases, is exposed to the 
surface. Drops form at nucleation sites which are probably 
randomly distributed. They are estimated to be on the order 
of one hundred million sites per square centimeter [2] . The 
nucleation size drops are of the order of one-tenth of a 
micron in diameter. These small drops initially are distri- 
buted over the entire surface and begin growing very 
rapidly due to condensation of vapor. When neighboring drops 
touch, they coalesce into a single new drop. If a nucleation 
site is exposed by the coalescence activity, a new nucleation 
size drop appears. The process of growth and coalescence 
continues until a drop reaches a size at which the surface 
tension forces holding the drop on the surface are exceeded 
by the gravity force on the drop mass. Drops of the depart- 
ing size are observed to slide down the surface and sweep 
other drops with them. This departing action also exposes 
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In 1930 Schmidt, Schurig, and Sellschopp [3] first 
investigated dropwise condensation as a means of heat transfer. 
Since their investigation, there has been a controversy as to 
the origin of the drops and as to whether the heat is 
transferred through the drops or through the bare surface 
between the drops. The opinions can be divided into two broad 
theories. First, the "film theory" contends that vapor 
condenses on the bare surface between the drops and that the 
resulting liquid forms into drops. Variations of this theory 
are supported by the work of early investigators such as 
Jacob [4], Eucken [5], and Emmons [6]. More recent investi- 
gators have proposed the "nucleation theory" which disagrees 
with the "film theory." The proponents of the "nucleation 
theory" contend that drops nucleate from discrete and 
randomly distributed sites on the condensing surface. 
Furthermore, the theory states that there is no liquid 
between the drops and that the path of heat transfer is 
through the drops only. In 1935 Tammann and Boehme [7] 
were the first of several to observe that drops form at the 
same sites on successive condensation cycles. In 1964 Umur 
and Griffith [8] proved that there is no liquid between the 
drops. Using thermodynamic considerations and an optical 
study, they demonstrated that a film no greater than a 
monolayer thick could exist between the drops on the surface. 
Umur and Griffith's work is final strong proof of the 
validity of the "nucleation theory." 
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In agreement with the nucleation theory, Rose 
[9,10] was the first to develop a method of predicting the 
heat transfer performance of dropwise condensation by 
analyzing the heat transferred through a single drop and then 
calculating the heat flux through the condensing surface by 
assuming a drop distribution. By fitting four constants into 
his equations , Rose was able to predict values which were 
comparable to experimental results. Continuing with a 
similar approach, Graheun [2] , further developed an analytical 
expression for the heat flux through a drop which included 
the heat transfer resistances of drop curvature, interfacial 
resistance between the liquid and vapor, and conduction. In 
order to authenticate his model, he experimentally measured 
the steady state drop distribution. The smallest drops 
which Graham measured and counted were ten microns in 
diameter. Drops of smaller diameter could not be accurately 
observed because of a resolution limitation of the optics. 

The small drops are moving very fast and have diameters on the 
order of the wave length of visible light. By combining his 
experimental distribution and his analytical model, Graham 
found that fifty percent or more of the heat flux was 
transferred through drops smaller than ten microns in 
diameter and concluded that his heat transfer model could 
not be verified unless the drop distribution was known 
over the entire range of drop sizes. Since a complete drop 
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distribution cannot be experimentally determined, Graham 
recommended computer simulation of the dropwise condensation 
process as a means of obtaining the complete drop distribution. 

The first computer simulation of steady state 
dropwise condensation was carried out by Gose, Mucciardi, 
and Baer [11] . The model included nucleation, drop growth due 
to condensation, coalescence, and removal of drops of the 
departing size. Their expression for drop growth from con- 
densation considered only the conduction resistance to the 
heat flux through the drop. For a departing drop with a 
diameter of two millimeters and for a site density of 
5 X lO'^ sites per square centimeter, their model predicted 
a heat transfer coefficient of approximately 460 British 
Thermal Units per foot square hour degree Fahrenheit.* This 
value is more than an order of magnitude lower than experi- 
mentally obtained coefficients. 

A later study by Tanasawa and Tachibana [12] used 
a similar technique and predicted a coefficient of 3,700 for 

3 

a site density of 4 x 10 . This prediction is lower than 

the value of 40,000 measured by Graham [2] at the same 

condensing conditions. Both of the previous computer 

simulations modeled drops of the entire range of sizes on 

a fixed area. This approach and their computer facilities 

*For this thesis heat transfer coefficients will have the 
dimensions of British Thermal Units per foot square hour 
degree Fahrenheit and site densities will be expressed as 
sites per square centimeter. 
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limited their investigations to artifically low site 

densities. Gose suggested that theoretical values would more 

closely match experimental results if larger site densities 
7 8 

of 10 to 10 could be modeled. These two computer simula- 
tions and Graham's work provided the essential background 
for this thesis. 

Using Graham's model of heat flux through a single 
drop, a computer simulation of the dropwise condensation 
process was carried out for this thesis in three parts. 

First, a computer model was constructed which predicts heat 
transfer coefficients and drop distributions for a given set 
of condensing conditions. Secondly, a parametric analysis 
was conducted varying the significant parameters, site 
density, saturation temperature of the vapor, and the tempera- 
ture difference between the vapor and the condensing surface. 
The third and final part of this thesis consisted of a compari- 
son of experimental results and model predictions. 
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CHAPTER II 
COMPUTER SIMULATION 



The Computer Model 

To model a representative condensing surface of one 
square centimeter would require the use of approximately one 
hundred million nucleation sites. The computer storage 
necessary to identify the location and activity of these 
nucleation sites was not available with the IBM 360-65 Model 
Computer with a Fortran G Compiler used for this study. If 
core storage were available, the computational time for high 
site densities would be large and the procedure would be 
uneconomical. To solve the storage constraint and to build 
an accurate and efficient computer simulation, an innovative 
modeling technique was developed to model the dropwise con- 
densation process on an initially drop free surface. 

The simulation models a single cycle in which a drop 
of the departing size is formed on an initially bare surface. 
To start the modeling process, the surface is exposed to an 
unlimited supply of vapor. Small primary drops appear at 
randomly distributed nucleation sites. These drops increase 
their size by condensation of vapor on their surface and by 
coalescence with neighboring drops. As the clock time 
increases, the average drop size existing on the surface 
increases. Eventually, a drop of the departing size is 
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produced and slides off the surface. The simulation process 
stops when the first drop of departing size is formed. 

The approach of the simulation is one of "viewing’’ 
the condensing surface with decreasing magnification as the 
average drop size on the surface increases. With a micro- 
scope one can see smaller objects by increasing the 
magnification, but by doing so, the field of vision is 
decreased. Similarly, the first area of interest is 
"viewed" with the highest magnification and is, therefore , the 
smallest area modeled. As the clock time increases, the 
average drop size increases. When the area of interest is 
no longer representative, the magnification is decreased and 
hence a larger area is visible. Ensuing drop growth 
necessitates lower and lower magnifications until departing 
drops are visible. An area of interest for a given magnifi- 
cation is called a stage. Table 1 shows how the stages are 
arranged for a representative calculation. Drops of the 
smallest size are not visible on later stages, but their 
presence is accounted for. 

The modeling technique starts with two assumptions. 
The simulation models the continuous process with discrete 
time intervals or "time steps". Drop growth from condensation 
and drop coalescence are modeled to occur sequentially within 
the time step. Whenever the simulation clock time is increased 
by the time step, a routine to determine drop growth is 
carried out, then the coalescence activity is modeled by a 
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separate routine. This procedure is considered to give an 
accurate representation of the actual continuous process. 
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TABLE 1 

ARRANGEMENT OF STAGES 



Stage 

Number Area in Square Microns 
1 10 ^ 

2 10 ^ 

3 10^ 

4 10^ 

5 10^ 

6 3.24 X 10® 



Range of Radii of "Visible" 
Drops in Microns 

0.10 to 3.75 

1.59 to 12.04 

8.03 to 40.64 

27.10 to 137.17 

91.45 to 462.96 

308.64 to 1250.00 



Note: This stage arrangement is for a saturation temperature 

g 

of 212 “F and for a nucleation site density of 10 
sites per square centimeter. 
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The growth routine determines a drop's size from 
an expression of drop age as a function of drop diameter. 

The age is increased by the time step. The drop's new size 
due to condensation during the time step is calculated using 
a relationship for drop diameter as a function of drop age. 

A detailed discussion of the growth routine is included in 
Appendix A. 

Upon completion of the growth routine, a check for 
coalescence is carried out. The coalescence routine employed 
by this study is innovative and consists of a method of 
selectively searching for coalescences between drops and for 
the covering and uncovering of nucleation sites. Each drop 
is considered, in turn, from left to right on the stage area as 
the focal point of the search. A search square of two drop 
diameters on a side is centered on the drop of interest. 

Only nucleation sites and drops with centers inside the 
search square are considered. A determination is made 
whether or not the central drop blankets any of the active 
sites or coalesces with any of its neighbors. If an 
uncovered site is found to be blanketed, it is made inactive. 
If the drop of interest is found to be touching a neighbor 
drop of the same size or smaller inside the square, they 
coalesce to form a single new drop at the central drop's 
location. After the coalescence occurs, a check for the 
uncovering of inactive sites is made. Sites uncovered by the 
coalescence are made active. 
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The coalescence has increased the size of the 
central drop. Consequently, an enlarged search area must be 
considered so that all possible coalescences of the central 
drop with neighboring drops of the same size or smaller are 
found. After each coalescence, the search area is enlarged. 
When the central drop fails to coalesce, the search recommen- 
ces around the next drop to the right on the stage area. 

When all drops have been the center of a search area, the 
coalescence routine ends. 

The central drop only coalesces with drops of the 
same size or smaller. By following this procedure for all 
search squares, all possible coalescences are found. This 
search technique has the advantages of requiring a small 
amount of core storage and is reasonably efficient for one 
thousand or fewer sites per stage area. Appendix B contains 
a more detailed explanation of the coalescence routine. 

To start the first stage requires selection of a 
time step value. If the time step is too large, the drops 
nucleated on the first time step will overlap by a large 
amount. Drops are observed to coalesce upon first touching. 
Significant overlapping is a serious departure from nature 
at high site densities. The time step for a given site 
density should be small enough to produce primary drops 
which do not artifically overlap. 
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The time step can also be too small. The problem 
arises because the approximation for drop radius as a function 
of age becomes inaccurate for very small drop radii. For this 
reason, the time step should be large enough to produce 
primary drops which are three to four times the minimum drop 
radius CAppendix A) . 

The First Stage 

O 

For a site density of 10 , the first stage has one 
thousand randomly distributed nucleation sites on a square of 
thirty-one microns on a side. On the first time step of the 
first stage, the growth routine causes hemispherical drops of 
uniform size to appear centered on nucleation sites. The size 
of the uniform drops is commensurate with the time step. 

Next, the coalescence routine is carried out. If neighboring 
drops touch, they coalesce to form a single new drop centered 
on the nucleation site of the formerly larger drop. At the 
end of the coalescence routine, the drops are no longer of 
uniform size. Their sizes vary depending on the number of 
coalescences which they have undergone. The stage surface 
contains nucleation sites which were uncovered by the 
coalescence activity as well as drops of different sizes. 

Other nucleation sites are inactive because they are 
covered by a drop. Upon completion of the coalescence 
routine, the simulation continues by increasing the clock 
time by the time step. 
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During the second time step, new drops nucleate 
and grow on the nucleation sites uncovered during the 
previous coalescence routine. Drops already existing on the 
surface grow by condensation. The amount of growth depends 
on the drop's size and the magnitude of the time step. The 
coalescence check is made, and the clock is again increased. 
The growth and coalescence sequence with increasing time 
continues until the stage termination criteria is reached. 

The stage termination criteria prevents "border 
effects". Drops near the border or edge of the stage area 
have an abnormally low nxamber of neighbors; therefore, they 
have less opportunity to increase their size by coalescence. 
The effect of this unrepresentative activity is minimized by 
modeling drops on a given stage only until they become 
relatively large when compared to the surface area. There 
are two tests in the program for stage termination. If a 
single drop is formed which covers greater than ten percent 
of the stage area, the stage ends. A category of smaller 
drops may also cause termination if they cover more than 
twenty percent of the stage area. These drops have radii 
which are forty percent smaller than the radius of the single 
drop which can cause termination. 

The amount of heat transferred through the 
condensing surface can be determined from the amount of 
liquid contained in the drops on the surface. The amount of 
heat transferred during a time step is simply the total 
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volume change of all drops on the surface multiplied by the 
density of the liquid and the latent heat of transformation. 
Since the rate of heat transfer, condensing surface area, 
and temperature difference between the vapor and the surface 
are known; the heat transfer coefficient for each time step 
can be calculated. From the total amount of liquid mass on 
the surface at a given clock time, the time averaged heat 
transfer coefficient for the process can be obtained and 
is the primary output of the stage. Figure 1 illustrates 
the time averaged heat transfer coefficient as a function of 
time. 

The Second Stage 

Since the area under study in the second stage is 
much larger than the first stage, the nucleation cannot be 
modeled in detail as it was on the first stage. Rather, it 
is assumed that on the bare area, nucleation and growth for 
the time interval equal to t^^ results in the same heat flux 
per unit area as growth of the first stage from time zero to 
time t^. Thus the heat flux for the bare area is calculated 
as the product of bare area, time interval t^^, and average 
heat transfer coefficient at time t^^; see Figure 1. This 
heat flux results in the creation of an amount of liquid on 
the bare area. Drops of radius r 2 are placed on uncovered 
locations on the bare area to account for this amount of 
liquid for each time step of the second stage. 
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The starting conditions of the second stage are 
determined by an interative matching procedure which matches 
the results of the first stage with the start up of the 
second stage. The start of the second stage is at the time 
t^ Csee Figure 1 ) which is half the cycle time of the first 
stage. The initial drops of the second stage are assumed to 
be uniform in size. A drop radius, r2» is assumed. The 
number of drops, N2, is calculated so that the total amount 
of liquid on the surface at time t^^ is the same for the first 
and the second stage; this makes the average heat transfer 
coefficient of the first and second stage equal. 

The N2 drops are distributed on a uniformly spaced 
net of locations, the number of locations, being larger 

than N2« At each location a random nvunber is used to deter- 
mine if a drop is present. All N2 drops are located on the 
locations. The locations are uniformly spaced to 
ensure that none of the drops coalesce immediately after 
they are placed on the surface. 

After all of the drops have been placed on the 
surface, the time is increased to 2 t^^ and the drops grow by 
condensation. Additional drops of radii r2 are placed on 
uncovered locations to account for the condensing productivity 
of the invisible drops on the apparently bare area. The 
average heat transfer coefficient at time 2 t^^ is calculated 
and compared to the value found from the first stage. If the 
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values disagree, a new radius X 2 is assumed and the process 
is repeated until the average heat transfer coefficient of 
the first and second stages agree at times t^^ and 2tj^. 

When the proper value of r 2 is obtained, the second 
stage is allowed to proceed with time steps equal to tj^. The 
process of growth and coalescence of existing drops and the 
creation of new drops on the bare area is continued for stage 
two until the termination criteria is met. At this time, the 
third begins with a larger area and uses a starting process 
which is analogous to the start up of the second stage. 

Time steps equal one half the cycle time of the second stage 
are used. CFigure 2 illustrates the start up of the third 
stage). During the growth of the third stage, bare area is 
assigned an average heat transfer coefficient found from the 
results of the second stage for the same time interval. 

Later Stages 

Following stages are all similar to the third 
stage. They use the same effective coefficient concept and 
termination criteria. Each stage differs from the former by 
having a lower magnification and having a stage area which 
is ten times the larger than the previous stage area. The 
simulation finally ends on the last stage. 

The final stage area is slightly greater than a 
factor of ten larger than the proceeding stage area. The 
size chosen is three and one quarter square centimeters and 
is comparable to the size of the experimental condensing 
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surfaces [2,13]. The termination of the last stage is 
different from the criteria used for the earlier stages . 

The last stage shuts off when a drop of the departing size 
is first formed. The values of the departing drop size used 
in this simulation were experimentally determined by Graham 
[2] . The final output of the complete simulation is a time 
averaged heat transfer coefficient for the process. The 
simulation program is visualized by a flow diagram. Figure 3. 
For the program, the input variables, source listing, and 
description of output are presented in Appendix C, Appendix 
D, and Appendix E, respectively. 
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TIME (SECONDS *10”^) 



FIGURE 1 FIRST STAGE OUTPUT 



THIRD STAGE DATA POINTS 
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FIGURE 2 THIRD STAGE START UP 
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FIGURE 3 



PROGRAM FLOW DIAGRAM 
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CHAPTER III 

DISCUSSION OF RESULTS 
Effect of Site Density 

At a saturation temperature of 212 ®F and at a 
temperature difference of one-half degree, heat transfer 
predictions were made for various site densities. 

Figure 4 displays the time average heat transfer coefficient 
over a complete simulation cycle for a range of site 
densities. The process is controlled by two limiting 
factors. Initially, the surface is covered by small drops 
which provide little resistance to the heat flux. They 
limit the process of heat transfer because they are in 
insufficient number to provide an adequate amount of drop 
surface on which vapor can condense. As time increases, the 
average drop size increases and the eimount of drop surface 
available for condensation increases. This effect accounts 
for the initial increase of the time averaged coefficient. 

The heat transfer coefficient reaches a maximum 
where another limiting factor becomes controlling. As the 
average drop size increases, the average drop resistance to 
heat transfer increases. After the maximum, the average 
drop resistance becomes the controlling limitation of the 
process. At higher site densities, the process is initially 
less limited by the inadequate amount of drop surface area 
and becomes resistance limited more rapidly. 

The higher site densities enable the process to 
obtain greater time averaged heat transfer coefficients. 
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Higher coefficient values are indicative of faster rates of 
mass accumulation. For this reason, drops of the departing 
size are produced sooner for process with higher site 
densities . 

Once the process becomes resistance limited, the 
rate of decline in the time average value of the heat trans- 
fer coefficient is the same for all site densities. This 
observation is based on the similarity of the slopes of the 
curves after the maximum has been reached. The decrease in 
the heat transfer coefficient after the maximum appears to 
be independent of site density. 

At the saturation temperature of 88 ®F and for the 
temperature difference of one-half degree, similar heat trans- 
fer predictions were made and are plotted on Figure 5 . At 
the lower saturation temperature, the interfacial resistance 
to heat transfer is greater by an order of magnitude and has 
the effect of decreasing the growth rate of smaller drops by 
a similar amount. The reduced growth rate at lower satura- 
tion temperatures makes the process more strongly limited 
by the initial inadequacy of available condensing surface. 
Consequently, at the lower saturation temperature smaller 
maximums occur later in time when the site density is held 
constant. 

A calculation of the maximiom theoretical heat 
transfer coefficient based on a maximum drop growth rate 
was made. The details of the calculation are given in 
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Appendix P. The values obtained are plotted on Figure 6 
which is a graph of the time averaged heat transfer 
coefficient at the end of the cycle for different site 
densities. The site densities examined produced coefficients 
well below the theoretical maximum. For a given site den- 
sity, a larger value of the coefficient was obtained at the 
higher saturation temperature because of the faster growth 
rates. As the site density increases, the difference in 
growth rates has a more pronounced effect and the coeffi- 
cient values for the two saturation temperatures diverge. 

As the site density decreases, the effect of the difference 
in growth rates diminishes and the coefficient values 
converge . 

Next, the time averaged distribution of drop 
sizes over the complete cycle for various site densities 
were plotted and compared. Figure 7 and Figure 8 show the 
number of drops per square centimeter for a given drop 
radius at the saturation temperature of 88®F and 212®F. At 
either saturation temperature and for drop radii greater than 
ten microns, the distributions are almost invariant with 
site density. For radii smaller than ten microns, the 
distributions diverge. At the higher site densities, there 
are more small drops. The presence of greater numbers of 
small drops account for the larger heat transfer coefficients 
predicted at higher site densities. 
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FIGURE 7 VARIATIONS OF THE DROP DISTRIBUTION 
WITH SITE DENSITY AT 88*F 
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DROP RADIUS (MICRONS) 



FIGURE 8 VARIATIONS OF THE DROP DISTRIBUTION 
WITH SITE DENSITY AT 212 ®F 
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Estimation of Experimental Site Densities 

A determination of actual site density was made by 
comparirg the model's drop distribution for the simulation 
cycle to Graham's measured distribution of the steady state 
process. Figures 9 and 10 display the model's prediction 
and Graham's data points for approximately the same heat 
transfer coefficient. The theoretical distributions agree 
well with Graham's measured distributions. Since the model's 
distributions are invariant with site density in the region 
measured by Graham, the actual site density had to be 
inferred by matching the model's predicted coefficient with 
Graham's measured heat transfer coefficient. For a tempera- 
ture difference of one-half degree, the model indicated that 
Graham's smooth copper condensing surface had an active site 

7 

density of about 2.5 x 10 sites per square centimeter at 
212®F and of about 2.0 x 10^ at 88®F. 

A standard experimental method of determining site 
density is to photograph a portion of the condensing surface 
with a large number of just visible drops. These drops have 
diameters close to the wave length of visible light; there- 
fore, the picture's quality is limited by the resolution of 
optics. By estimating the size of the drops and their num- 
ber, a drop density can be calculated. The nucleation site 
density can be assumed to be equal to or greater than the 
measured drop density. The model indicated that this 
technique may yield an incorrect estimate at high site 
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densities. For various nucleation site densities of the 
model. Table 2 shows the model's highest instantaneous 
drop density for different radii. Even if drops with a 
radius of one micron could be accurately measured, the 
nucleation site density could not with certainty be 
inferred. 
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FIGURE 9 COMPARISON OF PREDICTED AND 



MEASURED DROP DISTRIBUTIONS AT SS^F 
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DROP RADIUS (MICRONS) 

FIGURE 10 COMPARISON OF PREDICTED AND MEASURED 
DROP DISTRIBUTIONS AT 212°F 
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TABLE 2 

HIGHEST INSTANTANEOUS DROP DENSITY 



Model's Nucleation 

Site Density 1,0 

10® 8.1 X 10® 

lo"^ 3.6 X 10® 

10® 7.6 X 10® 

10® 1.0 X 10^ 



Drop Radius in Microns 

2.5 5.0 10.0 
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Note: The site densities are in dimensions of sites per 

square centimeter. The drop densities have units 
of drops per square centimeter. The predictions 
are for the condensing conditions of one-half 
degree temperature difference and for a saturation 
temperature of 212 ®F. 
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Effect of Temperature Differences 

At the saturation temperature of 88 °F and 212 °F and 
for the site densities determined by matching Graham's 
measured heat transfer coefficients, predictions were made 
for various temperature differences between the vapor and 
the surface. The results are displayed in Figure 11 and 
Figure 12, The model predicts an approximately constant 
heat transfer coefficient for a constant site density in the 
temperature difference range from a quarter degree to four 
degrees . The experimental results also plotted on the 
figures agree with the constant heat transfer coefficients 
values for temperature differences greater than one degree. 
Below one degree the heat transfer coefficients are 
experimentally observed to decrease. This disagreement 
between the model prediction for low temperature difference 
may be attributable to the model's assumption of constant 
site density. 

The minimum drop radius increases rapidly as the 
temperature difference is reduced. Figure 13 is a graph of 
the minimian drop radius plotted as a function of temperature 
difference. The curve indicates that as the temperature 
difference decreases from one degree to a quarter degree, 
the minimum drop radius increases by a factor of five for 
both 88°F and 212°F saturation temperatures. At high site 
densities , large values of the minimum radius could cause 
artificial overlapping on the first stage of the model and 



SATURATION TEMPERATURE = 88 ®F 
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FIGURE 12 VARIATIONS OF THE HEAT TRANSFER COEFFICIENT 
WITH TEMPERATURE DIFFERENCE AT 212®F 
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TEMPERATURE DIFFERENCE (®F) 

FIGURE 13 VARIATIONS OF MINIMUN DROP RADIUS 



WITH TEMPERATURE 
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thereby predict an erroneously high heat transfer coefficient 
The maximum theoretical site density for a given 
temperature difference was determined to see if the con- 
stant site density assumption at small temperature 
differences was causing artificial overlapping. Appendix G 
contains derivation of a maximum site density for a given 
temperature difference. Figure 14 is the result of a calcula 
tion using the derivation. Even for a temperature difference 
of a quarter degree, the theoretical maximum density is well 
above the site densities assumed for this analysis. The 
model's predictions at temperature differences less than one 
degree are believed to be valid for an assumption of constant 
site density. Kinematics of drop growth cannot account for 
the measured decrease in the heat transfer coefficient. 

The observed decrease can be explained by a reduction in the 
actual site density as the temperature difference decreases 
below one degree. 

Effect of Saturation Temperature 

To isolate the effect of saturation temperature, 
calculations were made for the saturation temperatures of 
88®F, 150®F, 212®F, and 274*F at the site density of 10^ and 
at the temperature difference of three degrees. Figure 15 
displays the theoretical prediction for the heat transfer 
coefficient as a function of saturation temperature at the 
constant site density. Above 150 *F the coefficient remains 
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FIGURE 14 MAXIMUN THEORETICAL SITE DENSITY AT 212®F 



TEMPERATURE DIFFERENCE = 3.0®F 
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FIGURE 16 VARIATIONS OF THE DROP DISTRIBUTION WITH 

SATURATION TEMPERATURE AT CONSTANT SITE DENSITY 
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approximately constant. Below 150 ®F the heat transfer 
coefficient decreases . 

The cause of the decrease can be understood by 
examining the drop distributions for the four saturation 
temperatures. Figure 16 displays the drop distributions. 

As the saturation temperature increases/ the number of 
small drops increases. This effect can be explained by 
considering the interfacial resistance. As the saturation 
temperature increases, the interfacial resistance decreases; 
and the drop growth rate of the smaller drops increases. 

For the constant site density assumption, it is the increased 
growth rate that is responsible for the greater nuitibers of 
small drops which appear in the drop distribution. Higher 
heat transfer coefficients are predicted when the drop 
distribution contains more small drops. 

Graham's experimental data for heat transfer 
coefficients as a function of saturation temperature is also 
plotted on Figure 15. The experimental result has a greater 
slope than the constant site density curve. This difference 
of slopes is further indication that the actual site density 
decreases as the saturation temperature decreases from 212 ®F 
to 88“F. 

Time Necessary to Produce a Drop of the Departing Size 

At a saturation temperature of 212 “F and for a 
temperature difference of two degrees, predictions were made 
by varying site density in order to match Westwater's data [13] 
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for heat transfer coefficients and sweeping frequency of 
drops of the maximvun size. The frequencies were converted 
into the time necessary to produce a drop of the departing 
size. Figure 17 shows measured and predicted heat transfer 
coefficients plotted as a function of the time necessary to 
produce departing drops. For the same heat transfer 
coefficient, the model predicts times which are twice as 
large as Westwater observed. Figure 18 shows the model's 
time averaged coefficient as a function of simulation time 
at the same condensity conditions. Reducing the model's 
cycle time by a factor of two means only a twenty percent 
increase in the model's predicted coefficient. For this 
reason, the lack of agreement with measured times is not 
considered to invalidate the other results of the model. 

This disagreement may be attributed to the 
hemispherical drop assxmiption or to the model's failure to 
account for sweeping. Drops of all sizes were assumed to be 
hemispherical in shape. A flatten or distorted drop of the 
same volume as a hemispherical drop will cover a larger 
portion of the surface. The larger surface coverage for the 
same drop mass implies a greater opportunity for coalescence. 
Since large drops grow mostly by coalescence [2,11], the 
distorted drops should mature faster than the perfectly 
hemispherical drops. 
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PRODUCE A DROP OF THE DEPARTING SIZE 
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The simulation is a model of the start of condensing 
on an initially bare surface and does not include the effect 
of sweeping. Experimentally, drops of the departing size 
slide down the surface and sweep vertical strips of the sur- 
face bare. In these bare regions, drops grow and appear 
approximately uniform in size. Since a drop can depart at 
any time in the actual steady state process, adjacent strips 
contain uniform drops of different generations. The regions 
of drops of different generations interact and may hasten the 
production of drops of the departing size. The model's 
failure to account for sweeping or the hemispherical 
assumption may be the reason for the overestimate of the 
time necessary to produce a drop of the departing size. 
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CHAPTER IV 
CONCLUSIONS 

1. Graham's model of heat flux through a single drop 
is a valid representation of the physical phenomenon. 

2. As the site density increases, the heat transfer 
coefficient increases because of a greater number of small 
drops in the distribution. 

3. For a temperature difference greater than one 
degree, the measured heat transfer coefficient remains 
constant because the site density is constant. For tempera- 
ture differences less than one degree, the decrease in 
measured coefficients can be explained by a reduction in the 
site density. 

4. As the saturation temperature decreases from 212®F 
to 88®F, the active nucleation site density decreases, 

5. The model predicts a time necessary to produce a 
drop of the departing size which is twice as large as is 
observed for the steady state process. 

6. This computer simulation provides a method of 
accurately predicting the heat transfer properties of the 
dropwise condensation process. 
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CHAPTER V 

SUGGESTIONS FOR FUTURE WORK 



1. The present model can be improved by using finer 
time steps and more elaborate starting technique for the 
later stages. Incorporation of the sweeping effect should 
also be added so that the simulation will more closely model 
the steady state process. 

2. The effect of non-hemispherical drops may require 
investigation. A combination of experimental observations 
and a theoretical analysis should determine the effect, 

if any, of non-hemispherical drops. 
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APPENDIX A 

MODEL OF DROP GROWTH DUE TO CONDENSATION 



Graham [2] proposed an analytical model for the 
growth of a single drop due to condensation of vapor. The 
heat flux associated with the condensing of vapor passes 
through three independent series resistances. The heat 
transfer resistances included are curvature, interfacial, 
and conduction. Each of these resistances is discussed 
below . 

1. Curvature Resistance. At a curved liquid-vapor 
interface, the saturation temperature of the vapor in 
equilibrium with the drop is less than the saturation 
temperature of the vapor in equilibrixim with liquid at a 
flat interface at the saime pressure. A derivation of the 
temperature drop through the interface due to curvature can 
be found in Graham's thesis. The result of his work is 
reproduced below: 



AT 



c 



2TgG 



1 

r 



(A.l) 



where AT^ temperature drop due to curvature 
Tg saturation temperature 
o surface tension 

latent heat of transformation 
p liquid density 
r drop radius 
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When the total temperature difference available between the 
vapor and the condensing surface is equal to the temperature 
drop due to curvature, the minimum drop size which can 
exist on the surface is determined. The expression for 
minimum drop radius is written below: 

2^s"^ 1 

^min HfgP AT^ 

where AT^ the total temperature difference between 
the vapor and the condensing surface. 

By combining the two previous equations, an expression for 
the temperature difference due to curvature can be found. 

ATc = (A. 3) 

where AT^ the temperature difference due to 
curvature . 

This temperature difference is not an actual resistance to 
heat transfer, but can be thought of as an equivalent 
resistance. 

2. Interfacial Mass Transfer Resistance. At the liquid- 
vapor interface, a pressure difference is necessary to force 
the mass transfer. The pressure difference can be converted 
into a temperature difference and is written below: 



AT^ = 



hj^2Trr 



(A. 4) 
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Where Q rate of heat transfer 



interfacial heat transfer coefficient 
AT^ temperature drop due interfacial 



resistance. 



The interfacial coefficient, hj^, can be calculated from the 
equation below which was derived by Nabavian and Bromley [14] . 



Figure A-1 is a plot of the interfacial heat transfer coeffi- 
cient as a function of saturation temperature. 

3. Conduction Resistance. The most significant 
resistance associated with the liquid drop is the resistance 
to heat conduction through the drop. The conduction resistance 
of 




T V 
s g 



(A. 5) 



where a condensation coefficient (assumed to be 
unity) . 



M molecular weight 
R universal gas constant 
Vg specific volume of the vapor 



R 



AT 

cyl Q 



KTrr 



r 



2 



(A. 6) 



INTERFACIAL HEAT TRANSFER COEFFICIENT, h. CBTU/HR*FT 



-64- 




FIGURE A-1 VARIATION OF INTERFACIAL 
HEAT TRANSFER COEFFICIENT 
WITH TEMPERATURE. 

REPRODUCED FROM GRAHAM [2] . 










% 






this shift. The age discrepancy is initially about 4 x 10 
seconds and decreases as the drop radius increases. The 
error induced by the imprecise knowledge of the nucleation 
size should not be significant because little vapor is 
condensed in the small time interval. 

Each resistance has been expressed as a temperature 
drop. The total temperature difference between the vapor and 
the condensing surface is equal to the s\am of the three 
resistance induced temperature drops. Then by utilizing the 
idea that the rate of heat transfer through a single drop 
is proportional to the drop volume increase with time , 

Graham obtained an expression for drop age as a function of 
drop diameter. 

Graham's derivation is reproduced below starting 
with the summation of temperature drops. AT^ is the total 
temperature difference between the vapor and the surface. 

AT^ — ATq + ATj_ + AT^c 

(curvature) (interfacial) (conduction) 

(A. 8) 

AT - -5^ AT + ^ (A. 9) 

^ ^ h^2Trr 4KTrr^ 

The heat flux through the drop is assumed to be proportional 
to the change in volume of the drop. 
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where Q rate of heat transfer 

AT temperature across the cylinder 

conduction resistance of the cylinder 
K thermal conductivity of liquid 
r drop radius 

For hemispherical drops, the above expression is multiplied 
by a shape factor of less than one. Mikic [15] determined 
the value of the shape factor to be one fourth. 

Multiplying by the shape factor and then solving 
for the temperature drop through the hemispherical drop, an 
expression for conduction resistance is obtained: 



AT 



dc 



Qr 

4Krrr^ 



(A. 7) 



where AT^^ temperature drop through the drop due 
to conduction resistance. 

Figures A-2 and A-3 are reproduced from Graham's thesis and 
indicated the effect of the three resistances on drop diameter 
as a function of drop age. 

Graham [2] , stated that nucleation radius was 
probably between one and two times the minimum drop radius. 

A drop of the minimum size has a theoretical zero growth 
rate. Several nucleation radii were examined to determine 
their effect. The curve of drop radius as a function of drop 
age kept essentially the same shape, but was shifted to the 
right as nucleation radius decreased. Figure A-4 illustrates 
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FIGURE A- 2 DROP GROWTH RATE FOR 88 *F. 

REPRODUCED FROM GRAHAM [ 2 ] . 
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FIGURE A- 3 DROP GROWTH RATE FOR 212 ®F 
REPRODUCED FROM GRAHAM [2] 
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FIGURE A- 4 VARIATIONS OF DROP GROWTH RATE 
WITH DROP NUCLEATION RADIUS 
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where V drop volxune 

t time or drop age 

dV/dt change in drop volume with respect to time 
For an hemispherical drop, 



dV 

at ~ 



2irr'‘ 



dt 



(A. 10) 



Q = PHfg 



Substituting equation (A. 11) 



AT. = at. + 

t r t h. 



into equation (A. 9), 



dr . ^^fg rdr 
dt 2K dt 



CA.ll) 



(A. 12) 



Solving for dr/dt and substituting drop dicimeter D for twice 
the drop radius , 



dD _ 4 

at PBfg 



AT. 



1 - ^min/° 
D/2K + 2/h. 



(A. 13) 



Equation (A. 13) can be integrated and then evaluated at the 
lower limit. Since the logarithmic terms become undefined at 
drop age zero when the drop diameter equals the minimxim 
diameter, the function must be evaluated at the minimum 
diameter plus epsilon. Graham's final result of drop age as 
a function of drop diameter is listed below. 
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>] 




1 




(A. 14) 



The simulation growth routine required an expression 
of drop radius as a function of age. For all cases in this 
study, Graham's equation (A, 14) was inverted by substituting 
a constant equal to five minimum drop diameter for the dia- 
meter variable in the logarithmic terms. This action con- 
verts the logarithmic expressions into constants. The 
inverstion is then finished by completing the square. The 
final result (A. 15) is an analytical approximation for drop 
diameter as a function of drop age and is shown below: 




(A. 15) 



where a = 1/4K 





2D„.„/h.) in (4D„.„) 



min' 



min 



d = 4AT^/pHfg 
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The approximation has reasonable accuracy over almost all of 
the range of drop sizes. The greatest inaccuracy experienced 
for the condensing conditions of this study was about two 
percent for drop radii of less than one micron. As the drop 
radius increased above one micron, the inaccuracy rapidly 
diminished. 

Using equations (A. 14) and (A.15) , the drop growth 
routine can be carried out. First a drop's age is determined 
using equation (A. 14). The age calculated is then increased 
by the time step. Finally, the new size of a drop is 
calculated with expression (A.15) . 
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APPENDIX B 
COALESCENCE ROUTINE 

The coalescence routine developed is capable of 
determining coalescence of drops, blanketing of drop sites 
or locations, and the uncovering of sites. Four variables 
associated with each site are required for the routine to 
operate with either randomly or uniformly distributed sites. 
The first two of the four site variables are coordinates, 

"X" and "y". These are obtained by superimposing a first 
quadrant coordinate system on the stage area. The coordin- 
ates of sites are necessary to calculate the distances be- 
tween sites. Then the sites are ordered from left to right 
on the stage square and assigned index numbers represented 
by the variable "I". The site with the lowest X coordinate 
value is site number one; the site with the highest X 
coordinate value has the largest index number. The fourth 
variable, ”R", identifies the activity at the site. The 
"R” variable can assume three types of possible values; 
minus one, zero, or a finite positive number. The minus 
one value indicates that the site is covered or blanketed by 
a drop centered on another site. The zero value represents 
a site that is exposed to the vapor or uncovered. The third 
"R" value is a finite positive number which is the radius of 
a drop in feet centered on the site. Touching of drops and 
covering of nucleation sites is determined by comparison of 
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the distance between sites and the drop radius. These 
four site variables; "X, Y, I, and R" ; are necessary to 
perform the bookkeeping required by the coalescence activity. 

The coalescence routine commences for each time 
step at the site with the lowest index niimber with a drop 
centered on it and proceeds until the site with a drop and 
the highest index number is examined for possible 
coalescences. A search area around the "Ith" drop of 
interest is defined as a square of two "Ith" drop diameters 
on a side. Figure B-1 is a diagram of a search square on the 
simulated condensing surface. Sites with a zero "R" value or 
a positive "R" value inside the square are examined. A 
determination is made whether or not the "Ith” drop with its 
increased radius from condensation growth blankets any un- 
covered or active sites. At the same time, a check for the 
possible coalescence between the central "Ith" drop and a 
neighboring "Jth" drop with a center inside the square is 
conducted. Only neighboring drops of the same size or 
smaller are not examined even though they may be blanketing 
part of the search or even touching the central drop. 

By maintaining a consistant examining procedure as 
the routine moves from drop to drop and from left to right on 
the simulated surface, all possible coalescences are found. 
This approach requires a relatively small amount of computer 
core storage and is reasonably efficient for one thousand 
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FIGURE B-1 SECTION OF STAGE AREA WITH 
SEARCH SQUARE CENTERED ON 
ITH DROP 
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sites on a stage area. For large numbers of small drops / 
small search areas are used. For the small number of large 
drops, large search areas are employed. This flexibility 
is the principle reason for the routine's efficiency. 

When a site is found to be covered, its "R" value 
is changed from zero to minus one. If the "Ith" drop is 
found to be touching a neighboring "Jth" drop inside the 
search square, a coalescence occurs. The masses of the 
"Jth" and "Ith" drops are combined to form a new drop at the 
central location of the "Ith" drop. A check for uncovering 
of sites, previously blanketed by the "Jth" drop, is then 
carried out in a similar manner as the coalescence was 
determined. When all sites previously covered by the "Jth" 
drop have been found and uncovered, the search for 
coalescence starts once again by defining a new larger 
search square around the "Ith" drop. The increased size of 
the square is indicative the the increase in size of the 
"Ith" drop due to the coalescence with the "Jth" drop. 

When no more coalescences are obtained between the "Ith" 
drop and any of its neighbors inside the search square, the 
coalescence routine centers its attention on the drop of 
the next higher index number. The same search for 
coalescence, covering, and uncovering is carried out. 

When the drop with the highest index number fails to touch a 
neighboring drop, the coalescence routine ends. 
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APPENDIX C 

DESCRIPTION OF THE INPUT DATA 

Seven data cards are used to read in fifteen 
variables necessary to run the program. These variables will 
be discussed in the order which they are read in. 

1. First Data Card 

Variable Name : IX 

Format; 19 

Description; The variable "IX" is any number with 
nine or less integers so long as the last integer is odd. 

The number is necessary to call a random number generator 
which is the IBM Scientific Subroutine RANDU. The random 
number generator is used to create the random site location 
for the first stage. If fewer than one hundred sites are 
used for the first stage, the selection of the variable "IX" 
can cause the time averaged heat transfer coefficient at the 
end of the first stage to vary by as much as ten percent. 

2. Second Data Card 

Variable Names: NN, SCALE 

Format: IlO, F10.4 

Description: These two variables determine the site 

density. "NN" is the number of sites for the first stage. 
"SCALE" is the dimension of the side of the first stage 
square. Scale is read in units of microns and then conver- 
ted into feet in the program. (The program does most of its 
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calculations in English Engineering Units.) Any number of 
"NN" sites can be used, but the scale can be only 10.0, 
31.623, 100.0, 316.23, or 1000.0 microns. Site densities up 
to ten to the ninth sites per square centimeter can be 
accommodated. (The program ran more efficiently by not 
letting the number of sites exceed one thousand on the 
initial stage area.) 

3. Third Data Card 

Variable Names: MM, TSTEP 

Format: 110, Ell. 4 

Description: "MM" is the number of time steps 

allowed for each stage. A limit of one thousand should not 
be exceeded. "TSTEP" is the time step for the first stage. 
The value selected should be large enough so that drops 
originating from nucleation sites have a radius which is 
three or more times larger than the minimum drop radius. The 
growth models inaccuracy for very small drops dictates this 
constraint. TSTEP has the units of seconds. 

4. Fourth Data Card 

Variable Names: TSAT, SURTE, HEATL, THCON, DTEMT, 

FDEN, HTCI 

Format: 6F10.4, Ell. 4 

Description : 

"TSAT" is the saturation temperature of the 
vapor for a flat vapor-liquid interface in degrees Rankine. 
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"SURTE" is the surface tension of the liquid 
at the saturation temperature and has the units of pounds 
force per foot. 

"HEATL” is the latent heat of transformation 
of the vapor to the liquid at the saturation temperature. 

Its units are British Themal Units per pound mass. 

"THCON" is the thermal conductivity of the 
liquid at the saturation temperature and has the units of 
British Thermal Units per foot-hour-degree absolute. 

''DTEMT" is the temperature difference between 
the vapor and the base of the drop which is approximated by 
the average surface temperature. The units are degrees 
Rankine . 

"FDEN" is the liquid's density in pounds mass 

per cubic foot. 

"HTCI" is the interfacial heat transfer 
coefficient and was obtained from a graph of values for 
various saturation temperature in Graham's thesis [2] . The 
coefficient has the units of British Thermal Units per 
square foot-hour-degree absolute. 

5. Fifth Data Card 

Variable Name ; RNUC 
Format: F6.3 

Description: "RNUC" is the size chosen for the 

nucleating drops necessary for the math model of condensation 
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growth. Graham [2] stated that the nucleation size was 
probably between one and two times the minimum diameter. 

This simulation used a nucleation drop radius of one and a 
half times the minimiam radius. The "RNUC" value for this 
case is the factor one point five. 

6. Sixth Data Card 

Variable Name: RMAX 

Format : FIO . 4 

Description: "RMAX" is the radius in microns of a 

drop of the departing size. 

The values used in this study were determined by 
Grahcun [2] . For a temperature difference of one-half 
degree, Graham measured values of 1500 microns and 1250 
microns for the saturation temperatures of 88 and 212 
degrees Fahrenheit respectively. Values of 1125 and 1375 
microns were assumed for the saturation temperatures of 
274 degrees and 150 degrees. (For this study, the maximum 
drop size was assximed to remain constant for all temperature 
differences . ) 
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APPENDIX D 

PROGRAM SOURCE LISTING 

The computer simulation of this study consists 
of a single program written in Fortran IV, Level G, for 
use on the IBM Model 360/65 Digital Computer. The follow- 
ing pages of this appendix contain a source listing of the 
program. 
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CALCULATION OF PA ARRAY 
PA( l) = TTIME(l»/TTIMEnR) 



- 101 - 



fv) 















A 






# 














m 






X 




















o 


































w 






QC 














UJ 


A 




UJ 














X 


oc 




a 
















»«4 


















»- 


•k 




uo 




















a 














m 


II 




o 














•> 






QC 














X 






O 














(\i 




















v-4 






U. 














•k 


1—4 




o 














m 


w 




















UJ 




ct 












«-4 


Hk< 


X 




UJ 












+ 








OQ 












CD 


< 






X 












»«4 


UJ 






r> 














QC 


o 




z 


A 










QC 


< 


•k 






ct 










o 


CD 


UJ 
















K 


CD 


«k 


— k 


< 












O 


f» 


W X 


QC 


h- 


UJ 










< 


•k 


< c 




3 


X 










u. 


X 


UJ «-4 


•k 


O 
















Qt: •- 




<I 


K 




or 








*-4 


< ^ 


II 




K 










»->4 


•k 


CD • 


• 4 


o 






•» 








• 


CO o 


««iik 


K 








> 




CD 


A 


•-« 


►-< «k 


X 






II 


< 






l»4 


► UJ 


^ — k 








*— • 


a: 




w 


W 


kwl » 


ct ^ 




1 










< 


< 


X 


O w 


> 


►— 




•k 


< 




UJ 


UJ 


^ C 


H cc 


< 


w 




— k 




fSI 


0^ 


QC 


< f-4 


U O 


QC 


UJ 


< 


»-p4 




ro 


< 


<. 


UJ »> 


< K 


QC 


X 


Ql 


w 


O 


'C 


CD 


< 


a: •4' 


u. u 


< 


»— « 




<r 


— ‘ K 




CD 


m 


< • 


• < 




K 


•» 


Q. 


'•t U 


c 




•k 


<r c 


^ IL 


i/) 


*- 


X 




• < 


K — 




X 


#-4 


X k-k 


t>0 


1 


»-4 


•k 


C IL 


CD 


CO 




•k UJ 


•-4 •» >0 


3 



ct w 



w LL IL 

o 



o 

o w 
< 



• w X 

o 



•»4 ^ • 

^ w C 
fvj 



> 
< 
C£, 
Ct 
< 

2T ^ 
ai 

UL O 



\rs 





UJ 






»-4 




X 


•k 




A 


Ok CL 




A 


1-4 




f -4 


4 «k »-4 


— k 


CL 


o 


2 CD 


•k 


X 




CO 


• 


r- 


O 2 


»-k 4 






II II 




h- 






•k 


^ - 


m 


•k 






•. o 


(VJ 


(—4 




CD 


•k 


m 


•• O 


II 










•—4 


•k 


o c 


f-* ^ *-4 


C 


2 


UJ 


^ • 


II 


K 




>o 


X 


vO 


kO •- 


»o4 




• 


CD *-4 CC 






X 


O »-i 


>o X 




f-4 


U 


M 


II o 


1—4 


K 


LU 


•k 


in 


•k 


»-4 ^ 




1 


K 


•-4 m »-4 


UJ 


•k 


o 


•k 


»-4 


•k OS 


•k 


1-4 


•—4 


»-4 


•-4 II UJ 




W 


D nO 


w 




w < 




•-4 


O 


W kO w 


r:> 


nO 


w 




w 


o ^ 


sO 


W 


LO 




— z> 


r: 


II 




w 


K 


•»— ' 


K ^ 




♦ 


• 


QC QC 


2 


w 


K 


k— r 


H 


w K 


w 


K 


QC 


< 


2 


m 




»->4 


UJ 


< 


UJ 


<X D 


CD QC 




O O m 




LU 


< 


UJ 


< 


UJ < 


LU 


< 


LU 




C^ w ^ 


<? 




K 


K 


X 


K 


X u 






kk -4 


H 1 - 1 - 


K 


►- 


X 


K 


X 


1 - X 


K 


X 


> 


K lA Z 1 - 




w 


2 


1 —* 


QC 


1-4 


QC ^ 




II 


W 


o u zr 


1—4 


QC 


1—4 


QC 


^ QC 


1—4 


QC 


2 




-J 2 


n < 


O 


QC 


C 


QC 


O < 


a 


CD 


CL 


< o < 


O 


oc 


O 


QC 


O 


QC O QC 


o 


O 


2 O UJ O 


O Ok 


O 


3 u. 


3 : LL O 


o 




»i >4 


U- O u. 


u 


:b 


Ok 


:s 


IL 


:k IL 3 r 


u. 


u 


»«4 


Q O U 






o 




GO 












Csl 


4-k4 




r^ 




K 


sO 




m 






»->4 






cn 




fD 




CD 








fO 


m 




*->4 




O 












C7> 






























o 












tr> 



o o 



u o o 



u o 



CALCULATION OF DELN ARRAY 
L0S=1 

KRV=KVV(L0S ) 

593 CONTINUE 

KRW=KKV(LOS) 
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APPENDIX E 

DESCRIPTION OF OUTPUT 

This appendix is added to assist a possible future 
user of the program. The output will be categorized and 
discussed in the order which it appears. 

1. The input variables defined in Appendix C are 
written as the first output. These units are given. Next, 
the result of a calculation of minimum drop radius and the 
variables for the drop growth model are written out. This 
information is contained in the first twelve lines of 
output . 

2. Each time step has ten lines of output. The time 
step output variables are written below as they appear. 

"HTC" is the instantaneous heat transfer 

2 

coefficient for the time step in BTU/HR FT ®F. 

"TIME" is the simulation clock time in seconds. 

"THTC" is the time averaged heat transfer 
coefficient. 

"COALESCENCE COUNTER" is the number of coalescences 
which occurred during the time step. 

Distribution categories are listed. "COVERED" 
implies that the drop site is blanketed or covered. 

"PRIMARY" implies that a drop site is uncovered. Other 
distribution categories are listed as drop radii in microns. 

A radius category is defined by a radius differential of 
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forty percent. Beneath the distribution categories, the 
numbers of each kind is listed. 

"TOTAL NO. OF DROPS" is simply the total nximber of 
drops on the stage surface at the end of the coalescence 
routine . 

"AVERAGE DROP SIZE" is the average size of a drop 
existing on the surface at the end of the coalescence routine. 

"DAREAM" is the amount of area covered by drops 
after coalescence and has the units of square microns . 

"PERC" is the percentage of the stage area covered 
by drops of the category which can cause termination of the 
stage. 

"BAREAM" is the amount of bare area in square 
microns after coalescence. This variable is the last of the 
time step variables. 

3. The third major kind of output data described is 
the start up procedure between the stages. The variables 
will be listed and described as they appear. 

"SCALE" is the dimension of the side of the next 
stage square in microns. 

"TSTEP" is the time step for the next stage in 

seconds . 

"THTCl and THTC2" are the time averaged heat 
transfer coefficients obtained from the pervious stage 
result. The start up iteration matches these values for the 
next stage. 
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Iteration sequence usually consist of several lines 
on which the iteration variables repeated. "LC” is the number 
of the iteration trial. "NRl" is the number of drops to be 
located on the stage area for the first time step. ”rmm 1'' is 
the radius of uniform drops produced on the bare area. "HCK” 
is a comparison of the time averaged heat transfer coeffi- 
cients for the two consecutive stages at the second time step 
of the later stage. "RMM12” is the radius in microns which 
the uniform drops obtain after growing for one time step. 

The iteration ends when the "HCK" value is within acceptable 
tolerances. The number of sites on the uniformly spaced 
stage area is indicated by the variable "NNLT". The 
appearance of this variable implies the completion of the 
start up process. 

4. At the end of the last stage, the output for the 
last time step contains the time averaged heat transfer 
coefficient for the entire process CTHTC) . The time to 
produce a drop of the departing size is listed as the 
variable "TIME". Also the size and number of departing drops 
are listed. 

5. The last significant output is an array containing 
the information of the drop distribution. The output 
variables are described as they appear. 

"WSSCD" values are the radius categories of the 



drop distribution. 
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"BWSS(I)" values are the time average distribution 
of drops resulting from the separate stage calculations. 

The units are drops per square centimeter. 

"DELNCD" values are the time average drop 
distribution for the complete simulation cycle and have units 
of number of drops per square centimeters. 

”HTCD(I}" is the ciamulative heat transfer 
coefficient for all drops of a given size ("WSS(I)”) or 
smaller. 

"THTCD" is the time averaged heat transfer 
coefficient for the entire cycle and is calculated from the 
distribution of drops. This output variable is the last of 
the program. 
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APPENDIX F 

MAXIMUM THEORETICAL HEAT TRANSFER COEFFICIENT 

An estimate of the maximum theoretical heat 
transfer coefficient can be made from the maximum rate of 
heat transfer through a single drop. Starting with an 
expression for the heat flux through a single, drop, the 
maximum theoretical coefficient is derived below. 



Q = pH. 



dV 



^fg dt 

For a hemispherical drop. 



dV _ 9__2 dr 

at - 



(F.l) 



(F.2) 



The heat transfer coefficient is defined as follows 



h - Q 1 
^ ~ A If 



(F.3) 



Substituting equation (F.3) into equation (F.4), 



h = 



2Trr2 

Wa 3F 



CF.4) 



The area of the base of the drop is irr' 
pH, 



h = 



fg 0^ 



Tf~ ^3t 



(F.5) 



Substituting in dD/dt for 2dr/dt, 

_ P^fg dD 
^ " ~If~ 3E 



(F.6) 



1 



r 
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The expression (F.6) for the heat transfer coefficient 
becomes a maximum when the drop growth rate variable 
(dD/dt) obtains its maximum value. Drop growth rate is 
defined by equation (A. 13) in Appendix A. Finally, the 
maximum theoretical heat transfer coefficient is expressed 
below as equation (F.7). 



h = /dD. 

max At '•3t ' max 



(F.7) 
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APPENDIX G 

DERIVATION OF MAXIMUM NUCLEATION SITE DENSITY ^ 

To ensure that the site densities used in this 
study were realistic, a derivation of a theoretical maximum 
density was made. It should be noted that the theoretical 
maximum is much greater than could be attained on a real 
condensity surface. 

Nucleation centers are assximed to be arranged in 
a triangular array, and the distance between centers is 
equal to twice the minimum drop radius (see Appendix A) . 

For a triangular array, ninety-one percent of the condensing 
surface is covered if the sites are circular. An equation 
can be derived which gives the maximum density as a function 
of the temperature difference and other properties. 

918 2 

Maximum Site Density = - — j ^ 10 (sites/cm ) 

irr 

(G.l) 

The radius, r, is equal to the minimum drop radius in 
microns and is defined below. 

_ 1 

'^min HjgP AT CG.2) 

Substituting in equation (G.2) into equation (G.l), 
the maximum theoretical site density as a function of the 
temperature difference and other properties can be obtained. 
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